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Point-contact (PC) investigations on the title compound in the normal and superconducting (SC) 
state (T c ~ 10.6 K) are presented. The temperature dependence of the SC gap of TmNi2B2C de- 
termined from Andreev- reflect ion (AR) spectra using the standard single-gap approximation (SGA) 
deviates from the BCS behavior in displaying a maximum at about T c /2. A refined analysis within 
the two-gap approximation provides evidence for the presence of a second gap twice as large as the 
main gap (the first one), while the latter is close to that within the SGA. This way, TmNi2B2C 
expands the number of nickel borocarbide superconductors which exhibit a clear multiband charac- 
ter. Additionally, for the first time "reentrant" features were found in the AR spectra for some PCs 
measured in a magnetic field. The PC spectroscopy of the electron-boson interaction in TmNi2B2C 
in the normal state reveals a pronounced phonon maximum at 9.5 meV and a more smeared one 
around 15meV, while at higher energies the PC spectra are almost featureless. Additionally, the 
most intense peak slightly above 3meV observed in the PC spectra of TmNi2B2C is presumably 
caused by crystalline-electric-field (CEF) excitations. The peak near 1 meV detected for some PC 
spectra is connected with a modification of the CEF probably due to boron or carbon vacancies, 
allowing to probe the local stoichiometry by PC spectroscopy. 

PACS numbers: 73.40.Jn, 74.45. +c, 74.70.Dd 



I. INTRODUCTION 



Compared with the isomorphic nonmagnetic transition 
metal borocarbides i?Ni2B2C with i?=Y, Lu, Sc, La, Th, 
the rare-earth nickel borocarbide -RM2B2C (here R is a 
rare-earth element) compounds with i?=Dy, Ho, Er, and 
Tm are of special interest due to two coexisting and com- 
peting ordered states in these systems: antiferromagnetic 
(AFM) and superconducting (SC), observed by electri- 
cal resistivity, specific heat, magnetization, neutron, etc. 
measurements (see Refs.— ~— and Refs. therein). The en- 
ergy scales of these two states are of similar size, with the 
SC critical temperature T c both larger (for R= Ho, Er, 
and Tm) and lower (for i?=Dy) than the AFM Neel tem- 
perature Tn- The SC state here is most probably of the 
BCS type and mediated by the electron-phonon inter- 
action. The magnetism in these rare-earth compounds 
is provided by their localized 4/ rare-earth ion's mo- 
ments which order at relatively high temperatures com- 
pared with their dipole-dipole interaction due to the indi- 
rect Ruderman-Kittel-Katsuya-Yosida (RKKY) interac- 
tion mediated by the itinerant electrons. Moreover, the 
crystal electric field (CEF) splits the degenerate 4/ rare- 
earth ion states. Thus, CEF effects may influence the 
formation of the AFM and SC states. For a deepened 
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understanding of .RN12B2C compounds spectral informa- 
tion both on the electronic structure and on the bosonic 
collective excitations are of importance. 

Point-contact (PC) spectroscopy^ is a unique tool to 
study simultaneously both the SC gap or order parameter 
and to provide spectral data as to normal state bosonic 
excitations (phonons, magnons, CEF excitations etc.) 
always present in these solids. Measuring the second 
derivative of the I(V) curves for the PC's, information 
about the electron-phonon interaction function a 2 F(u>)^ 
and analogously for other bosonic collective interactions 
in the normal state can be obtained. Also the behavior 
of the SC gap (i. e. of the order parameter in a more 
rigorous sense) in the SC state can be determined from 
the dV/dI(V) curves using a now routine approach 7 -. 

In this paper we focus on PC measurements on a 
rather special member of the rare-earth nickel borocar- 
bide family, TU1N12B2C, which becomes superconduct- 
ing below 11 K. But only below Tn — 1.5 K does an 
AFM ordering appear where transversely polarized spin- 
density waves with an incommensurate modulation of the 
magnetic moments sets inji— Thus TmNi2B2C exhibits 
the largest difference between T c and Tn in the whole 
.RNI2B2C series, which points to a relatively weak mag- 
netism which cannot destroy the superconductivity on 
large Fermi surface sheets (FSS's), but only weakens it, 
although sizably, strongly dependent on the purity of the 
single crystals and as well as on the field direction.— ~— 
The superconductivity is nevertheless affected by the spe- 
cial incommensurate magnetic state and its fluctuation 
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above Tn — 1.5 K, at least for somehow disordered sam- 
ples (probably with respect to C-B local disorder). In 
this context the short periodicity of the incommensu- 
rate antiferromagnetic structure, sw 25 A, as compared 
with the superconducting coherence length, £ « 125 A, 
is of interest^ and a further argument given there that 
magnetism and superconductivity affect each other. In 
the present paper we do not consider the very coexis- 
tence of magnetism and superconductivity because this 
is already outside the range of our experimental setup 
(T > 1.6 K). For the temperature range above T N and 
at least for low external magnetic fields one may also 
expect to detect remnants of the multiple-gap features 
well-known for non-magnetic borocarbides^— A clear 
hint for such a multiple-gap behavior is given by experi- 
mentally observed positive curvature pronounced here for 
the upper critical field applied || to the a6-plane. Due to 
the large magnetization for H || to the c-axis, H C 2 is most 
strongly suppressed and its positive curvature as wellJ^ 
To the best of our knowledge, no PC data concerning 
TmNi2B2C, are available, except our recent preliminary 
results! 14 ' 15 On the contrary, there are several papers de- 
voted to investigation of other members of the nickel 
borocarbide family (see Refsi 1 ^— and Rcfs. therein), 
where the SC gap as well as the electron-phonon spectral 
function were studied. The peculiar non-monotonic tem- 
perature dependence of the upper critical field H C 2(T) \\ 
to the c-axis observed frequently^ has been attributed 
to Pauli paramagnetic effects^. All this makes investi- 
gations of the missing TmNi2B2C by PC spectroscopy 
desirable in order to get a complete picture. Here we 
restrict ourselves to a study of the paramagnetic phase 
above Tn — 1.5 K. 



II. EXPERIMENTAL DETAILS 

TmNi2B2C single crystals were grown by the float- 
zone technique with optical heatingi 1 ^ For homogeniza- 
tion and in order to minimize the amount of possible 
disorder on the boron and carbon sites, the crystals were 
wrapped into tantalum foil and subjected to a heat treat- 
ment at 1000° C for 72 h and subsequently at 500°C for 
72 h under purified argon atmosphere. The residual re- 
sistance ratio [RRR =p(300K) / p(llK)] of the investi- 
gated crystal significantly increased by this heat treat- 
ment: from RRR=14 before heating to RRR=28 after it. 
The high perfection of the single crystals is confirmed by 
the small width (about 0.45 deg) of neutron diffraction 
peaks obtained for (103) Bragg reflection. The investi- 
gated single crystal had a resistivity of p (UK) = 1.7 
pfl cm in the normal state. Its SC transition tempera- 
ture amounts T c = 10.6 K, determined from the midpoint 
value of the SC transition. 

The PC spectra have been measured mainly along the 
c direction by a standard " needle-anvil" method^ touch- 
ing a cleaved TmNi2B2C surface with a sharpened thin 
Cu wire. A number of PC spectra has been measured 
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FIG. 1: (Color online) Normalized to the normal state and 
symmetrized dV/dl curves of a TmNi2B2C-Cu contact (R = 
5.7 $1) established approximately in the c direction at tem- 
peratures from the bottom to the top: 1.6 K, 2K, 2.5 K ... 
and further with 0.5 K step until 9.5 K (upper curve). Left 
inset: Example of a raw dV/dl curve at 1.6 K (solid curve) 
and 10 K (dashed curve) extended to the higher biases. Right 
inset: Some selected normalized and symmetrized dV/dl for 
the same PC (points) at T=1.6, 3.5 and 6K (from the bottom 
to the top) along with calculated theoretical (according to the 
theory^) curves (thin lines). 

perpendicular to the c direction, but we have found no 
qualitative difference between them and the PC spec- 
tra along the c direction. The reason might be that the 
cleaved surface was quite irregular and rough, therefore 
the crystallographic direction in PC is not well defined 
and the PC spectra are expected to be averaged over a 
large angle range. 

The differential resistance dV/dI(V) and the second 
derivatives d 2 V/dI 2 iV) were recorded through sweeping 
the dc current / on which a small ac current i was su- 
perimposed. The alternating voltage V\ oc dV/dI(V) 
as well as its second harmonic V 2 oc d 2 V/dI 2 (V) have 
been recorded using the standard phase sensitive lock- 
in technique. It is well known that the first deriva- 
tive dV/dI(V) = R(V) of the I(V) curve between 
a normal metal and a superconductor, the so-called 
Andreev-reflection (AR) spectrum, reflects the SC gap.— 
The second derivative of the /(V) curve of the type 
R~~ 1 dR/dV = R~ 2 d 2 V/dI 2 is directly proportional to 
the electron-boson spectral function according to the PC 
spectroscopy theory (see Ref£, chapter 3 and Ref£). 
This derivative can be readily obtained: simply by divid- 
ing the measured V 2 by Vf, i.e., R- l dR/dV cx V 2 /Vf. 

III. PC SPECTROSCOPY OF THE SC ENERGY 

GAP 

To obtain a spectroscopic insight into the SC energy 
gap as well as into the spectral function of the electron- 
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boson interaction, the contact size d should be smaller 
than the inelastic electron mean-free path. To avoid 
a variation of the SC gap on the scale of the PC size, 
also the contact diameter d should be less than the co- 
herence length £, which is in TmNi2B2C between 11 
and 15 nm.— i 17 i 20 The elastic electronic mean-free path in 
TmNi2B2C is estimated as about I ~ 20 nm, using typi- 
cal for the nickel borocarbides pi ~ 3.6 x 10 _12 f2 cm 2 (see 
Refj2!) and our residual resistivity po ~ 1.7/zilcm. Ac- 
cording to Wexler's formula— R PC — IGpl/S^d 2 + p/2d, 
(here we neglected contribution to the second term from 
the resistivity of a normal metal (needle) and suppose 
geometrically equal parts occupied by each metal in the 
case of heterocontact), the PC diameter d can be esti- 
mated from their resistance Rpc and pi. Using a typical 
Rpc between 1 f2 and 1017, one estimates a d value be- 
tween 30 nm and 10 nm. Thus, the necessary condition of 
the smallness of the PC size compared to I and £ may be 
fulfilled for the investigated PCs. Even more importantly 
from the experimental or practical point of view is the ob- 
servation of characteristic AR features in the dV/dI(V) 
spectra in the SC state along with the phonon structure 
in d 2 V/dI 2 (V) in the normal state, which can be used as 
criteria for a spectroscopic (ballistic or diffusive^) regime 
in PC. 

We succeeded in obtaining high-quality dV/dl PC 
characteristics as shown in Fig.[TJ which demonstrates 
clear AR features - deep minima around V — ±A/c 
at T <C T c , here A is the SC gap, and the absence of 
any other irregularities like spikes, humps, etc. For most 
PCs, the structure (minimum) which characterizes the 
SC state disappears at about 10 K, slightly below the 
mentioned bulk T c ~ 10.6 ± 0.05 K. The experimental 
dV/dl curves were fitted according to the well-known 
modified BTK theory (see Ref£, section 3.7) within the 
single gap approximation (SGA). In the inset of Fig.[TJ as 
an example, the calculated dV/dl along with the mea- 
sured ones are shown at several temperatures. From the 
fitting of experimental dV/dl at different temperatures, 
the temperature dependencies of the SC gap A, the so- 
called Dynes or broadening parameter the barrier 
strength Z, and the scaling factor S have been deter- 
mined (see Fig. [5] and Appendix). 

As shown in Fig. [21 T and Z are nearly constant, ex- 
cluding the region close to T c . However, in this region 
above 6 K the double minimum structure in dV/dl is 
washed out and the self-action computer procedure gives 
more room for the fitting parameters by approaching T c . 
The most interesting feature is the temperature depen- 
dence of the gap A(T). As shown in Fig. [21 A(T) devi- 
ates from the BCS-type behavior, reaching a maximum 
not at T -> 0, but at about T ~ 5K ~ T c /2. The major- 
ity of investigated PCs reveals a similar behavior for the 
extracted A. Moreover, the minima positions in dV/dl 
which at low temperatures should be close to the A value 
also display maximum (see Fig. [5]) directly supporting the 
obtained specific A(T) behavior in this temperature re- 
gion. The upper critical field H C 2(T) measured for fields 
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FIG. 2: (Color online) Upper panel: Temperature dependen- 
cies of the fitting parameters (within the SGA) for the PC 
spectra from Fig.[T] the SC gap A (triangles), the broadening 
parameter V (squares), the barrier parameter Z (diamonds) 
as well as the scaling factor S in the inset. The minimum 
position in dV/dl is shown by circles. Bottom panel: Phase 
diagram for both superconductivity and AFM spin-density 
wave order of TmNUI^C using data frorn^ along with our 
critical magnetic field data along the c-axis (squares) and ab- 
plane (circles). 



along the c direction and within the ab plane behaves 
similarly (see Fig. [2]), that is, it has a maximum. This 
broad maximum is closely related to the reentrant be- 
havior found in the temperature dependence of both the 
bulk resistivity (see inset of Fig. [5]) and of the magnetic 
susceptibility of TmNi2B2C in magnetic field. Obviously, 
the recovery of the normal state at low temperature in 
the presence of a sufficiently high applied magnetic field 
is responsible for the anomalous decrease of H C 2(T) at 
low temperatures. This H C 2(T) anomaly might support 
our finding. In fact, since in the clean limit one has 
in the orbital approximation (ignoring a possible Pauli 
limiting behavior mentioned abovei^) H C 2(T) oc A 2 (T), 
the anomalous decrease of the upper critical field is more 
pronounced in the gap itself. 

The gap distribution for the selected 14 contacts with 
pronounced AR minima along with the characteristic ra- 
tio 2A/keT c are shown in Fig. [31 The SC gap value Aq (in 



FIG. 3: (Color online) Distribution of the single SC gap A(0) 
extrapolated by the BCS curve to T=0 K (left panel) and the 
characteristic ratio 2A(0)/kBT c (right panel) for the measured 
contacts. Notice the different shapes of both histograms re- 
flect the different T c values of various probed PCs. 



the SGA extrapolated by the BCS dependence to T=0) 
results in a somewhat wide range Ao = 1.1 — 1.7 meV, 
which reveals a characteristic ratio 2A/kBT c between 
3 and 4.2^ here T c is the temperature where the AR 
features for the corresponding contact vanish, in other 
words, the local T c . It varies between 9 and 10 K. Then 
the mean 2Ao/ksT c exceeds the BCS weak coupling 
value 3.52, which points quite generally to a strongly 
coupled SC state. Alternatively, it might also reflect 
simply the local gap anisotropy on a given FSS. In or- 
der to explain the nearly isotropic H C 2 in DyNi2B2 C 24 i 25 
(see, e.g., Fig. 6 in Ref.—), a somewhat anisotropic 
electron-phonon interaction and also a corresponding gap 
anisotropy must be assumed to compensate the differ- 
ent anisotropy of the Fermi velocities of the calculated 
cushion FSS, which would also affect H C 2. A similar en- 
hanced gap ratio is also characteristic for other i?Ni2B2C 
compound s] 14 ' 15 measured by AR spectroscopy. The 
gap distribution in TmNi2B2C might be connected both 
with the anisotropy and the multiband superconductiv- 
ity, as suggested by the analysis of the PC data for the 
nonmagnetic LuNiaBzC^i 2 ! YNizBaCjii and the AFM 
ErNi2B2C— superconductors, or for the well-known two- 
gap superconductor MgB2<22 On the other hand, exist- 
ing tunneling spectroscopy data for TmNi2B2G2£ read: 
Ao = 1.45 meV and 2Ao/kBT c — 3.1 with error bars of 
± 15%. Additionally, they did not show a remarkable 
anisotropy and exhibit a BCS-like weak-coupling tem- 
perature dependence of the gap in the whole temperature 
range, even in the AFM region below 1.2 K. 

Analyzing in more details our data and the SGA fit- 
ting parameters, it turned out that: 1) usually a quite 
large T > A/ 2 is needed to fit the dV/dl curves, 2) the 
scaling factor S (see Appendix) varies considerably (see, 
e.g., Fig.[2J) , 3) there is a wide distribution of the gaps. 
The large T might point to a spatial (or directional) gap 
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FIG. 4: (Color online) Temperature dependencies of the fit- 
ting parameters in the two-gap approximation (TGA): small 
(closed triangles) and large (open triangles) the SC gaps Ai : 2, 
the broadening parameters ri = T2 — 0.2 (squares), the bar- 
rier parameter Z (diamonds), the scaling factor S (stars), and 
the contribution (weight factor) w to the dV/dl curve of the 
large gap (open circles). Error bars for the gaps show the 
difference between the presented gap values and the gaps cal- 
culated supposing T2 — 2Fi (see also footnote Ref. #54). 



distribution due to anisotropy or the multiband character 
of the SC state, a picture which is now widely accepted 
in the borocarbide community. Therefore, supposing a 
multiband SC state we have also carried out a two-gap 
fitj 5 ^ In general, this fit is similar to the one-gap fit shown 
in Fig. Upright inset), only the two-gap fit is a bit better 
in the 4-6 mV range for the dV/dl curves at the lowest 
temperatures. The results are shown in Fig.|4j First of 
all we note that 1) T is constant and more than two times 
smaller than that of our SGA fit 5 -^ and 2) the scaling fac- 
tor is close to 1 and varies moderately, except the region 
close to T c , where the fit gives more freedom for fitting 
parameters due to a temperature broadening and smear- 
ing of a double-minimum structure of dV/dl. These two 
observations suggest that the obtained results of our two- 
gap fit are physically reasonable. As it is seen from Fig. [4] 
the small gap exhibits a more or less similar behavior as 
the gap obtained in the SGA fit, while the second gap 
is two times larger, but its contribution to the spectra 
is about 10%, once again excluding the region close to 
T c . Thus we found plausible arguments for the presence 
of a second SC gap in TmNi2B2C. To clarify its nature, 
detailed directional PC investigations are desirable. 

In the end of this section we would like to mention that 
except YbNi2B2C with heavy-fermion like properties, at 
first glance the electronic structure as predicted by the 
first electronic structure calculations based on the local 
density approximation (LDA) is hardly affected by the 
rare-earth component directly^ In particular, all these 
compounds exhibit a special FSS, called in the case of 
i?=Lu, Y, Ho and Dy as "cushion"— or "pillow"^ which 
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is formed only by Ni 3d states, i.e., without an admix- 
ture of rare-earth 5d derived states. The latter medi- 
ate the RKKY exchange coupling of the magnetic mo- 
ments of the 4/ electrons to some of the conduction 
electrons whose superconductivity is then strongly af- 
fected by the presence of the localized rare-earth mag- 
netic moments. On the contrary, the superconductiv- 
ity on the "cushion" FSS is almost perfectly protected 
against local magnetic fields caused by the exchange in- 
teraction between its conduction electrons and the 4/- 
moments. However, recent band structure calculations 
for TmNi2B2C show explicitely that some FSS has some- 
what changed their shape. In particular, the "cushion" 
has not been detected.— By symmetry an analogous Ni- 
3d derived FSS free of an admixture with Tm 5d states 
should nevertheless exist. Since its exact shape has not 
been resolved yet, we denote it, for the sake of simplicity, 
as "pseudo-cushion" (pc)-FSS. On these grounds a sim- 
ilar local SC gap of ~ 1 meV can be expected for that 
pc-FSS in TmNi 2 B 2 C, too. The PC data for the smaller 
gap reported above (which is less sensitive to the vicinity 
of magnetism) provide strong support for such an as- 
signment in the title compound, too. The larger second 
gap reported here, too, with a stronger sensitivity to the 
vicinity of magnetism corresponds to the major FSS in 
nonmagnetic borocarbides, which in the present case is 
somewhat affected by the fluctuations of this AFM state 
occurring below 1.6 K. 



A. "Reentrant" features for dV/dl characteristics 
in an external magnetic field 

Let us turn to the measurements of dV/dl in an ex- 
ternal magnetic field. By studying the obtained curves 
one can distinguish two types of dV/dl behavior. One 
of them is shown in the upper panel of Fig. [5] Here the 
external magnetic field does not substantially change the 
position of the dV/dl minima (reflecting the SC gap) up 
to fields of about 1 T. The magnetic field reduces the in- 
tensity of the whole dV/dl minimum and its influence 
looks like a smearing of the AR minima. A qualitatively 
similar behavior is observed for most of PCs based on 
the title compound and on many other superconductors, 
including nickel borocarbide o 11 ' 12 ' 16 ' 35 as well as simple 
SC metals, e.g., Nb^ 6 .. In the meantime the dV/dl for 
some studied PCs have demonstrated a quite different 
behavior jS 5 - similar to that shown in Fig.(5Jb). Here, the 
double-minimum structure undergoes significant modifi- 
cations in a magnetic field. Below 0.6 T the dV/dl curves 
behave like in the previous case, but at higher field they 
transform into another structure with a plateau-like max- 
imum at zero bias, so that the zero-bias resistance reaches 
its normal-state value at 0.7 T and then this normal-state 
region spreads out to higher voltages by further field in- 
crease until all dV/dl curves transform into a parabolic- 
like normal-state behavior. 

This structure of dV/dl in Fig. [51(b) looks similar to 
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FIG. 5: (Color online) a) Raw dV/dl curves of contact from 
Fig.Q] (after the temperature measurements) measured at 
T=3K in magnetic field, b) Raw dV/dl curves of another 
TmNi2B2C-Cu contact measured at T=2K in a magnetic 
field. Note the transformation of the zero-bias maximum into 
a plateau- like structure above 0.7 T. Inset: reentrant features 
in the resistance of a TmNi2B2C single crystal in an external 
magnetic field B || c. 



reentrant behavior found in the temperature dependence 
of the bulk resistivity of TmNi2B2C in the external mag- 
netic fieldS 6 -, where the normal state is restored at 2 K in 
a field above 0.7 T (see inset in Fig.JSb). This behavior 
of the bulk resistivity, i.e., the recovering of the normal 
state at low temperature under a magnetic field, may be 
connected with some kind of magnetic ordering induced 
by the external magnetic field. The temperature increase 
(e. g., above 2K at 7T, see inset in Fig.[Sb) destroys this 
state and recovers the SC state. We assume that in our 
case the applied voltage plays a role similar to the tem- 
perature, namely, a relaxation of strongly nonequilibrium 
electron distribution function in a PC under an applied 
voltage 6 - can produce nonequilibrium phonons with an ef- 
fective temperature <eV/4 3 -Z. Therefore in our case ap- 
plied voltage increase recovers the SC state (SC features 
in dV/dl) in the field above 7T at 2K. 
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FIG. 6: (Color online) Magnetic field dependencies of the 
fitting parameters within the SGA taken at T—3K: the SC 
gap A (triangles), the broadening parameter F (squares), the 
barrier parameter Z (diamonds), the scaling factor S (stars), 
as well as the minimum position in dV/dl (circles) extracted 
from the PC spectra in Fig.[5ja). 
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FIG. 7: (Color online) R- 1 dR/dV oc d 2 V/dI 2 curves mea- 
sured at 2 K in a magnetic field for the 3.8 Q. PC from 
Fig-EJb). The inset shows the PC spectrum of another 22 Q. 
PC with sharper phonon maxima along with the phonon DOS 
for LUN12B2C (symbols). 



This qualitatively different behavior observed in both 
types of dV/dl characteristics is consistent with the ob- 
servation that reentrant behavior of p(T) in a magnetic 
field is not found for all TmNi 2 B 2 C single crystals^ 

We have also recovered the magnetic field dependence 
of the SC gap and of the other fitting parameters again in 
the SGA. The results are shown in Fig. [5] using the spec- 
tra from Fig.[5{a). It turned out that the effective single 
gap is almost magnetic field independent, strictly speak- 
ing, it exhibits even a slight increase. This observation 
is in accord with the minima position in dV/dl. How- 
ever, this unusual, absolutely counterintuitive and strik- 
ing behavior (never observed for any superconductor to 
the best of our knowledge) might be ascribed to an un- 
resolved anisotropic multigap scenario. In fact, provided 
the PC reflects also some (a6)-geometry, a gap contribu- 
tion from non-cushion FSS's with a much higher upper 
critical field i/ c2 than the depicted highest field of 1 T 
might readily resolve the puzzle. 

r steeply increases in an external magnetic field 
as it has been found for many other type-II 
superconductors^^ while in type-I superconductors T 
is almost field independent^ The increase of T is usu- 
ally connected with a shortening of the bosonic (Cooper 
pairs) life-time due to the pair-breaking effect of an exter- 
nal magnetic field4£ The scaling parameter S decreases 
with the magnetic field. This can be understood as a 
shrinkage of the SC volume in the PC region, probably 
due to vortex penetration. However, the determination 
of the SC gap and other parameters from AR spectra 
in a magnetic field using the standard BTK approach 
with the introduction of the SC density of states smear- 
ing due to life-time effects^ only is, in principle, a sim- 
plified approach. As mentioned in Refi 4 ^, the SC density 



of states in the mixed state varies in space, i. e., becomes 
nonuniform. This special situation requires an appro- 
priate theoretical description (not available at present) 
for the density of states to compute realistic AR spectra. 
Moreover, the local SC density of states in PC will de- 
pend on the position of the pinned vortex with respect 
to the PC core. Therefore the obtained effective param- 
eters of the homogeneous model and their magnetic field 
dependence extracted from AR spectra measured in an 
external magnetic field above the lower critical field H c \ 
should be interpreted with some caution. 



IV. PC SPECTROSCOPY OF BOSONIC 
EXCITATIONS IN THE NORMAL STATE 

Now we consider the measured PC electron bosonic 
interaction spectra (R~ 1 dR/dV = R~ 2 d 2 V/dI 2 curves) 
of TmNi2B2C-Cu PCs. These spectra show clear max- 
ima at about 3mV and 9.5 mV (see Fig. [7]) and a more 
smeared maximum at 15mV (see inset in Fig. [7J and the 
following figures). The maxima at about 9.5 mV and 
15 mV were not affected by magnetic field and corre- 
spond well to the low-energy maxima in the neutron 
phonon density of states of the nonmagnetic isostruc- 
tural compound LuNi 2 B 2 C 4 ^ (see Fig. [71 inset), only the 
high-energy part of the obtained PC spectra remains 
almost featureless. The displayed in the PC spectra 
phonon features are similar to those observed in other 
nickel borocarbide compoundsj 1 ^ 14 ' 16 ^ 7 ' 43 ^ 4 Thus, the 
observed peaks reflect the TmNi 2 B 2 C phonons excited 
by the inelastic scattering of electrons (energized by eV- 
bias). Unfortunately, even for the PC with a relatively 
large resistance (see Fig.[7J inset) for which the condi- 
tions for the ballistic (spectroscopic) regime is fulfilled 
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FIG. 8: (Color online) Upper panel: Position of the phonon 
peak on R~ 1 dR/dV oc d 2 V/dI 2 PC spectrum (shown in in- 
set) for several TmNi2B2C-Cu contacts versus magnetic field. 
Bottom panel: position of CEF peak(s) from Fig.[7](triangles) 
and Fig. [9] (circles). Dashed lines demonstrate linear fit of the 
data. Inset shows the typical PC spectrum with mentioned 
peaks marked by the arrows. 



with a larger probability, so far we could not resolve the 
phonon features above 20 mV, contrary to the case of 
HoNi2B2CJ^ One of the reasons might be that the PC 
spectrum of TmN^I^C contains a dominant 3-mV peak, 
which is connected with CEF excitations as is shown 
below. Thus, excitations of Tm ions to the first CEF 
level by e^-energized electrons followed by subsequent 
scattering on 9.5-mV and 15-mV phonons will lead to a 
shortening of the inelastic electron mean free path and 
to deviations from the ballistic (spectroscopic) regime in 
our PC with further increasing of bias eV. This blocks 
energy-resolved spectroscopy at higher energies. 

Next we will consider the mentioned low-energy peak 
near 3 mV, where the peak positions are marked by the 
tilted dashed line in Fig.[7J It can be assigned with the 
first excited CEF transition in TmNi2P>2C according to 
Ref.—. As shown in Fig.[H a magnetic field shifts this 
maximum to higher energies whereas the maximum at 
about 9.5 mV remains fixed. This confirms phononic ori- 
gin of the 9.5-mV maximum. The other CEF transitions 
at the higher energies mentioned in Ref.— could not be 
resolved in the PC spectra, probably due to their lower 
intensity and the mentioned deviation from the spectro- 



scopic regime with the bias increase. Thus, this low-lying 
first CEF excitation seems to play a dominant role in the 
transport properties of Tmf^P^C. The presence of such 
low-energy excitations in TmJ^P^C can be anticipated 
from the nonsaturated resistivity in this compound by 
lowering temperature with pronounced slope of p(T) un- 
til T c (see Fig.fj2», inset). 

We have also occasionally measured PC spectra with 
an additional peak at low energy near 1 mV (see Fig. [9]). 
A magnetic field shifts this peak as well as the second 
peak at about 3mV (see Fig. bottom panel) to the 
higher energies, whereas the positions of the phonon max- 
ima near 10 and 17 mV are stable. The 3-mV peak looks 
like the CEF peak on the previous spectra. The 1-mV 
peak persists up to the highest fields available in our 
experiment (7.5 T) exceeding the H c2 in the ab and c di- 
rections by a factor of four or two, respectively. Such 
a behavior in a magnetic field excludes any residual su- 
perconductivity in our PCs as the origin of this peak^. 
The appearance of this 1 mV peak may testify to the 
nonstoichiometry of TmN^P^C in the PC core. In con- 
trast to i?=Ho, Er compounds, here in the region of C 
and/or B vacancies or disordered C and B sites the orig- 
inal magnetic Tm-ion ground state doublet splits into 
two singlet states (see Ref4, page 768). If so, then on the 
one hand the PC spectroscopy enables us to distinguish 
the "fine structures" in the CEF schema of TmNi2B 2 C, 
while on the other hand it allows us to control the local 
stoichiometry. Another possibility is that TmNi2B2C ex- 
hibits a quadrupolar ordered phase below Tq=13.5K^. 
The consequence of this would be that the quadrupo- 
lar ordering lowers the CEF symmetry, resulting in the 
mentioned ground-state splitting. In this context we note 
that a weak peak at about of 1 meV has been observed 
in the inelastic neutron spectra of Tm .o5Yo.95Nl2B 2 C4£ 
It was ascribed by the authors to a transition derived 
from a splitted, originally two-fold degenerate Tm-ion 
ground-state. Thus, we also assume here a similar tran- 
sition resulting from a modified CEF Tm-ion ground 
state, either due to nonstoichiometry (vacancies) or due 
to a quadrupolar ordered phase. To finally elucidate the 
nature of these peaks, more detailed directional mea- 
surements in higher magnetic field with samples of well- 
defined stoichiometry are required as well. 



V. CONCLUSION 

We have carried out investigations of the SC gap and of 
the electron-bosonic interaction spectra in TmNi2B2C by 
means of AR and PC spectroscopy, respectively. The SC 
gap in TmNi2B2C is found to deviate considerably from 
the usual BCS-type behavior in showing a broad maxi- 
mum at T* ~ T c /2. An improved two-band analysis of 
our data also shows the presence of a second SC gap being 
about two times larger than the former, which remains 
close to the small "first" gap observed. The smaller gap 
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FIG. 9: (Color online) PC spectrum R~ 1 dR/dV of TmNi2B2C-Cu PC for two bias polarity with additional lmV peak vs 
magnetic field. Left inset: dV/dl for the PC from the main panel at low fields, which also demonstrates "reentrant" features 
above 500 mT. Right inset: PC spectrum for another PC at two fields showing a similar low-energy 1-mV peak structure. 



has been assigned to a FSS which has a similar orbital 
structure and a gap value as the "cushion" FSS found 
in other nickel-based 1221 borocarbidc compounds. This 
not-yet-resolved FSS for TmN^E^C is almost protected 
from the exchange interaction with the magnetic mo- 
ments of the rare-earth ions. The observed second larger 
gap is only a bit smaller than the larger gap in nonmag- 
netic borocarbide superconductors. This way TmNi2B2C 
interpolates between " cushion" -FSS-dominatcd almost 
single-band superconductors DyNi2B2C and H0N12B2C 
within the commensurate AFM state and the nonmag- 
netic borocarbide superconductors YN12B2C, LUM2B2C, 
and probably, also the SCM2B2C compound which, how- 
ever, is less well studied. 

We succeeded to measure PC electron-bosonic interac- 
tion spectra of TmNi2B2C with distinct phonon features 
at 9.5 mV and 15meV along with intense low-energy 
maxima of a nonphonon nature at about 3 meV and ad- 
ditionally at about 1 meV for some spectra. The 3-meV 
maximum is due to excitation of Tm ions on the upper 
(first) CEF level by eV-energized electrons. As to the 
1-meV peak, it is likely connected with a modification of 
the CEF due to boron or carbon vacancies. This allows, 
in principle, for the PC spectroscopy to be regarded as 
a local probe to check the stoichiometry. Intense peaks 
of CEF excitations in PC spectra point to their signifi- 
cant contribution to the spectral function of the electron- 
boson interaction. 
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Appendix 

As a rule, the parameter T is used to describe the 
nonthermal smearing of the experimental dV/dl curves, 
namely, the broadening of the characteristic AR minima 
in dV/dl, while the original physical meaning of T in 
the Dynes model^ as a finite lifetime of charge carriers 
due to inelastic scattering is usually not under consid- 
eration. This is because other effects contribute to the 
dV/dl broadening, e. g., as SC gaps distributio n 48 ' 49 (via 
nonhomogencity, surface, multiband, etc. effects) or gap 
anisotropy (see also §4.3.4 in Ref<^). 

Including T in the fit procedure leads not only to the 
smearing of the whole calculated dV/dl curve but also 
results in a reduced intensity of the whole AR (minima) 
structure. In this context we note, that for a high-quality 
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fit of dV/dl in full, not only the shape of dV/dl must be 
described, but also its absolute intensity must be correct. 
In the present paper we also address this point, namely, 
we considered the, so-called, scaling factor S, which is 
included to fit the intensity of calculated and experimen- 
tal curves. For instance, 5"=1 means that the calculated 
curve fits the measured dV/dl also in absolute values. 
In the case of the ordinary superconductor with a large 
coherence length like Zn, 39 we succeeded in performing 
such a fit practically for all PCs. If S" <1, the mea- 
sured dV/dl has a reduced intensity due to suppression 
of the AR signal or the opening of some non-AR chan- 
nel in the PC conductivity. Thus, the authors of Ref.— 
assumed that S <1 is realized when the Fermi surface 
is only partially gapped in the SC state. Of course, the 
more common explanation of non-SC regions in the PC 
or multicontact scenario, including PCs with suppressed 
superconductivity, cannot be ruled out as well. Thus, in 
Ref.— a contribution to the conductance from the nor- 
mal vortex cores in the superconductor was taken into 
account to describe a progressive suppression of the AR 
features in an external magnetic field. 



It turned out that even a seemingly nonphysical situa- 
tion with S >1 can also take place. As a rule, large V val- 
ues result in large S > 1 . A more reasonable explanation 
would be the following: if the smearing of dV/dl is due 
to a variation (distribution) of SC gaps over the Fermi 
surface, then the description of this smearing by includ- 
ing r results in a reduced calculated dV/dl intensity and 
finally leads to a scaling parameter S > 1. Hence such a 
situation with S > 1 might point to a gap distribution or 
to a multiband scenario. 

We should also mention that a contribution to the 
dV/dl minima may come from the Maxwell resistance^ 
which cannot be ignored for materials with a high resid- 
ual resistivity like, e. g., the heavy- fermion compounds*^ 
Since the residual resistivity in i?Ni2B2C compounds is 
above l/xf2cm, then the Maxwell contribution to Rpc can 
exceed 1S1 for a PC with a diameter of 10 nm. The van- 
ishing of the Maxwell resistance in the SC state can lead 
to a situation where the PC resistance (the conductance) 
can decrease (increase) more than 2 times, thus it exceeds 
the limit of the maximal AR signal. 
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